PCT/AU2004/000312 





PRIORITY 
DOCUMENT 

SUBMITTED OR TRANSMITTED IN 
COMPLIANCE WITH ROLE 17.1(a) OR(b) 



Patent Office 
Canberra 



^ k** t tj atypt? FXAMINATION SUPPORT AND 
I, JULIE BILLINGSLEY, ^.^SeProvisional specification 
SALES hereby certify that annexed -^^^q^901237 for a patent by THE 




WITNESS my hand this 
Thirtieth day of March 2004 



JULIE BILLINGSLEY 
TEAMT.RADF p t?V AMTNf ATION 
RTTPPORT A^Ti SALES 



BEST AVAILABLE COPY 



AUSTRALIA 
Patents Act 1990 

PROVISIONAL SPECIFICATION 

for the invention entitled: 

"Fractal Image Data and Image Generator" 



The invention is described in the following statement: 



P^OPEKUtABVaaa frwul Image dau prov.Ax-13/W/W 



-1- 

FRACTAL IMAGE DATA AND IMAGE GENERATOR 
FIELD OF THE INVENTION 

The present invention relates to fractals, and in particular both to a new class of fractals 
and aiso to a fractal generator able to iteratively generate image data and a number of such 
5 fractals. 

BACKGROUND 

Since the 1980's, fractals have been used to describe a wide variety of physical systems. 
Fractals can be generally classed as deterministic, if they can be generated by a 
deterministic process, or random, if the generation process uses one or more probability 

10 distributions to select iterated function systems. One of the best known applications of 
fractals is the simulation of natural formations occurring in the physical world, such as 
mountains, clouds, and waves. One of the difficulties of existing fractals and procedures 
for their generation is that deterministic fractals tend to be too rigidly deterministic, and 
random fractals too random, when these are compared with natural processes and forms 

15 occurring in the physical world. Deterministic fractals, while initially eye-catching, tend to 
lose the viewer's interest because they appear too repetitive. On the other hand, random 
fractals can also appear somewhat boring because the eye finds it hard to discern patterns, 
similar to the way that wispy white clouds can appear boring. These weaknesses are one of 
the factors that have limited the use of fractals in the content industries, in the areas of 

20 special effects in digital animation, computer games, and advertising, for example. 

It is desired to provide at least a useful alternative, and in particular to provide image data 
representing at least a nascent fractal and a generator and process for generating the image 
data that alleviate one or more of the above difficulties. 
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SUMMARY OF THE INVENTION 

In accordance with the present invention, there is provided a V-variable fractal 

5 In particular, the present invention provides a V-variable fractal represented by fractal 
image data, where V is an integer greater than one and represents the number of 
constituent images available for iterative combination to generate the fractal. The images 
are combined in a random manner. 

10 The present invention also provides image data, wherein a variable number n of constituent 
images are iteratively transformed and combined in a random manner to generate the 
image data, with 1 < n < V. 

Preferably, the constituent images are transformed using projective transformations at each 
15 iteration. 

The present invention also provides a fractal generation process, including: 

(i) randomly selecting images from a set of input images; 

(ii) selecting transformation functions from a set of transformation functions; 

20 (iii) generating transformed images by applying the selected transformation 

functions to the selected images; and 

(iv) generating an output image by combining the transformed images. 

Preferably, the process includes: 
25 (v) repeating steps (i) to (iv) to generate a set of output images; and 

(vi) repeating steps (i) to (v) using said set of output images as said set of input 
images to generate a new set of output images. 

Preferably, the process includes repeating step (vi) until said new set of output images is 
30 substantially independent of the first set of input images used in the process. 
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Preferably, the number of selected transformation functions is less than the number of 
transformation functions in said set of transformation functions. 

Preferably, the step of selecting transformation functions includes selecting an iterated 
5 function system from a set of iterated function systems, each iterated fiinction system 
including a set of transformation functions. 

Preferably, the selection of an iterated function system is based on selection probabilities 
associated with said iterated function systems. 

10 

Advantageously, the combining of said transformed images may include superimposing 
said transformed images. 

Advantageously, said transformation functions may include geometrical transformations. 

15 

Preferably, said geometrical transformations include scaling and translation. 

Advantageously, said geometrical transformations may include scaling, translation and 
geometrical distortion. 

20 

Preferably, said transformation functions are contractive on average. 

Advantageously, said geometrical transformations may be contractive transformations. 

25 Advantageously, said transformation functions may include projective transformations. 

Advantageously, said transformation functions may include transformations of brightness 
and/or colour. 

30 Advantageously, each of said transformation functions may be represented by one or more 
parameters. 
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Advantageously, the process may include generating said transformation functions and 
said selection probabilities. 

Advantageously, said transformation functions and said selection probabilities may be 
generated on the basis of predetermined probability distributions. 

The present invention also provides a fractal generation process, including randomly 
selecting from a set of input images, transforming the selected images, and combining the 
transformed images to generate a set of output images, and iterative repetition of these 
steps using the set of output images of each iteration as the set of input images for the next 
iteration. 

Advantageously, said selecting may include selecting the same input image more than 
once. 

Preferably, said transforming includes scaling and translating the selected images. 

Advantageously, said transforming may also include geometrically distorting the selected 
images. 

Preferably, the transforming is contractive. 

The present invention also provides a fractal generator, including: 

an image selector for selecting M of V input images; 

a function selector for selecting a set of M transformation functions; 

at least one image transformer for respectively applying the selected transformation 
functions to the selected input images; and 

a compositor for composing an output image from the images output by said at 
least one image transformer. 
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Advantageously, the fractal generator may include a set of V image buffers storing the 
input images. 

Preferably said function selector selects said set from N sets of transformation functions. 

5 

The present invention also provides a system having components for executing the steps of 
any one of the above processes. 

The present invention also provides program code for executing the steps of any one of the 
1 0 above processes. 

The present invention also provides a computer readable storage medium having stored 
thereon program code for executing the steps of any one of the above processes. 

15 The present invention also provides fractal data generated by any one of the above 
processes. 

The said fractal data represents one or more V-variable fractals. 

20 The present invention also provides a fractal generation system, including an image 
selector for selecting images from a set of input images, and an image transformer for 
transforming the selected images to generate a set of output images, said system being 
adapted to provide said set of output images as the set of input images to iteratively 
generate fractal image data. 

25 

Preferably, said image transformer includes one or more image transformation modules for 
transforming said selected images, and an image combination module for combining the 
transformed images. 

30 Preferably, said one or more image transformation modules are adapted to scale and 
translate said selected images. 
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Advantageously, said one or more image transfonnation modules may also be adapted to 
geometrically distort said selected images, 

5 Preferably, the system includes a transformation selection module to select transformations 
to be applied to the selected images from a set of transformations. 

Preferably, said transformation selection module is adapted to select transformations based 
on selection probabilities associated with said set of transformations. 

10 

The present invention also provides fractal image data representing a combination of two 
or more constituent first images, each of said first images representing a random 
transformed combination of two or more constituent second images, each of said second 
images representing a random transformed combination of two or more constituent third 
15 images, each of said third images representing a random transformed combination of two 
or more constituent fourth images, wherein each transformation includes at least one of 
translation and rotation. 

Preferably, each transformation is a projective transformation, such as an affine 
20 transformation. 

Preferably, each transformation includes contractive scaling. 
Preferably, the transformations are contractive on average. 

25 

Advantageously, each combination may be a superposition. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments of the present invention are hereinafter described, by way of 
example only, with reference to the accompanying drawings, wherein: 

Figure 1 is a block diagram of a preferred embodiment of a fractal generation 
5 system; 

Figure 2 is a block diagram of a fractal generator of the fractal generation system; 
Figure 3 is a block diagram of an image transformer of the fractal generator; 
Figure 4 is a flow diagram of a fractal generation process executed by the fractal 
generation system; 

10 Figure 5 is a schematic diagram illustrating transformation of an input image by a 

transformation module of the image transformer; 

Figures 6 to 10 are screenshots of images generated by the system using a first set 
of input parameters; 

Figures 11 to 15 are screenshots of images generated by the system using a second 
1 5 set of input parameters; and 

Figures 16 to 23 are screenshots of images generated by the system using a third set 
of input parameters. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

20 As shown in Figure 1, a fractal generation system includes a fractal generator 100, As 
shown in Figure 2, the fractal generator 100 includes an image selector 202, an iterated 
function system (BFS) selector 204, and an image transformer 206. The fractal generation 
system executes a fractal generation process that generates output image data 102 from 
input image data 104, iterated function systems 106, and selection probabilities 108. The 

25 input image data 104 constitutes a set of input images 104, and the output image data 102 
constitutes two or more output images 102 that may be considered as representing two or 
more fractals, as described below. 

In the described embodiment, the fractal generation system is a standard computer system 
30 such as an Intel™ x86-based personal computer including a Pentium™ processor 110, 
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random access memory (RAM) 112, and non-volatile (e.g., magnetic disk) storage 114, 
and the fractal generator is implemented by software modules stored on the non-volatile 
storage 114 of the computer and executed by the processor 110, However, it will be 
apparent that at least parts of the fractal generator can be alternatively implemented by 
5 dedicated hardware components, such as application-specific integrated circuits (ASICs). 

The fractal generation system is described as generating output image data 102 from input 
image data 104. Although it is not necessary that the input and output data relate to images 
as such, Le. 9 the system processes these as arbitrary numeric data, it is expected that the 
10 system will be used to process input data and generate output data that in each case will be 
visualised as one or more images on a display or hardcopy. Accordingly, the input image 
data 104 is hereinafter referred to as input images 104, and the output image data 102 is 
hereinafter referred to as output images 102. 

15 A user of the fractal generation system provides to the system a number V > 1 (i.e, at least 
two) input images 104, together with N > 1 iterated function systems F = {F ! , F N } 106 
and N selection probabilities P = {P ! , . . P N } 108. Each of the N iterated function systems 
106 comprises a set of M transformation functions. Each transformation function defines a 
contractive affine transformation that can be applied to an image, and is represented by the 

20 six parameters r, s, x, y, 0 and <p. Other transformations can be used, each represented by a 
set of parameters. As shown in Figure 5, an input image 502 (in this case having 
horizontal and vertical dimensions of unity) is transformed by first scaling the horizontal 
and vertical dimensions of the input image 502 by the factors r and s, respectively, with -1 
< r < 1 and -1 < s <1, followed by translation along the x and y axes by the amounts 

25 specified by the x and y parameters, respectively, followed by a geometrical distortion 
whereby horizontal lines are rotated by the angle 0, and vertical lines are rotated by the 
angle <p 9 with -180 £ 0 £ 180 and -180 <> <p <> 180 . In Figure 5, the angles 0 and <p, are 
nearly equal, and the output images 504 appears to retain its rectangular shape. However, 
when these two parameters are substantially unequal, a rectangular input image 502 is 

30 distorted to a generally parallelogram or diamond-like shape. 
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Because each transformation function is represented by the six transformation parameters 
described above, each iterated function system (IFS) 106, being a set of M transformation 
functions, can be represented as a parameter table 506 with M rows, each row representing 
a distinct transformation function. For example, an IFS 106 for M=3 can be represented as 
5 a 3x6 matrix or table, as follows: 



ri 


Si 




yi 


e, 


<pi 




S2 


x 2 


yi 


e 2 


<P2 


T3 


S3 


X 3 


y* 


e 3 


93 



Thus the iterated function systems 106 are provided as a set of N> 1 parameter tables, 
each having M > 1 rows. The N selection probabilities 108 are used to select a particular 
10 IFS parameter table) to apply to selected images, as described below. 

As shown in Figure 2, the fractal generator 100 operates on a set of V input buffers 208 to 
212, and a set of V output buffers 214 to 218 allocated from the computer RAM 112. 
Although Figure 2 shows V=3 of each kind of buffer, there can be in general an arbitrary 

15 number V > 1 of each buffer type, subject only to system constraints. The number of input 
images 104 provided to the system is less than or equal to the number V of input and 
output buffers, as the same image can be allocated to more than one buffer, and the input 
images 104 are copied into the respective input buffers 208 to 212. The fractal generator 
100 then executes a fractal generation process, as shown in Figure 4. The fractal generator 

20 100 generates output images into each of the output buffers 214 to 218 in a sequential 
manner, starting with the first output buffer 214, and ending with the last output 
buffer 218. When the outermost loop of the fractal generation process is executed, an 
iteration of the process is said to have completed, and occurs when all the output buffers 
214 to 218 have images stored in them. 

25 

The fractal generation process begins at step 402, when the first output buffer 214 is 
selected to receive the first output image. At step 404, the image selector 202 randomly 
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10 



selects M of the V input images stored in the input buffers 208 to 212, with the same image 
allowed to be selected more than once. Accordingly, each time a selection is made, the 
probability of selecting the input image from a particular buffer is equal to 1/V , although it 
will be apparent that alternative selection methods can be used. The M selected images are 
provided to the image transformer 206. At step 406, the IFS selector 204 randomly selects 
one of the N IFSs 106 using the respective N selection probabilities 108. For example, the 
N selection probabilities 108 can be provided as respective weights used to divide the 
numeric interval between 0 and 1 into N contiguous subdivisions, where the size of each 
subdivision is proportional to the weight for that subdivision. A random number between 
0 and 1 can then be generated to select one of the N IFSs 106, according to which 
subdivision the random number falls within. 



As shown in Figure 3, the image transformer 206 includes M transformation modules 302, 
304, and a superimpose module 306, shown for the example of M=2. At step 408, the 

15 transformation modules 302, 304 transform the selected input images using the selected 
IFS. The selected input images are transformed using respective transformation functions 
from the selected IFS. That is, the first transformation module 302 applies a first 
transformation function defined by the first row 312 of the selected parameter table 316 to 
the first selected input image 308. Similarly, the second transformation module 304 applies 

20 a second transformation function defined by the second row 314 of the selected parameter 
table 316 to the second selected input image 310. 

Each of the transformation modules 302, 304 generates a transformed version of its input 
image, using the corresponding transformation function, as described above. It will be 

25 appreciated that, although distinct transformation modules 302, 304 are shown in Figure 3, 
the M selected input images could alternatively be processed sequentially by a single 
transformation module. However, it may be particularly advantageous to provide distinct 
execution instances of each of the transformation modules 302, 304 if the system includes 
more than one processor 1 10 or if the transformation modules 302, 304 are implemented as 

30 dedicated hardware components. 
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In the preferred embodiment, the input images 104 provided to the system by a user are 
mxn rectangular images utilising a rectangular co-ordinate system with origin 0,0 
corresponding to one of the comers of the rectangular input images 104. The transformed 
images generated by the transformation modules 302, 304 are maintained as m x n images 
corresponding to the same spatial co-ordinates as the originai input images 104. 
Accordingly, a given transformed image will generally include a number of transformed 
versions of previous images positioned within an otherwise empty m x n space within the 
particular transformed image. Any transformed images or parts of transformed images that 
fall outside the m x n co-ordinate space of a transformed image are cropped and are not 
included within the transformed image. Because each transformed image corresponds to 
the same mxn region of coordinate space, transformed images can be combined by 
superposition. However, it will be apparent that alternative embodiments can be devised 
such that a transformed image is not generated as an m x n rectangular image within the 
mxn coordinate space defined by input images 104, and the combining of images may 
involve only partial overlap of the constituent images, or may not involve any overlap at 
all 

At step 410, the resulting transformed images are processed by the superimpose module 
306. This involves superimposing or overlaying the transformed images to generate a 
single output image 318. The superimposition is achieved by assigning to each pixel of the 
output image 318 a colour value equal to the average of the colour values of the 
corresponding pixels of the transformed images. The resulting superimposed output image 
318 is stored in the current output buffer, in this case the first output buffer 214. 

It will be apparent that the transformed images can be combined in a variety of alternative 
ways. For example, the images may be combined by putting one down on top of the 
preceding one, one after another, so that only the last value written to each pixel is kept. 
Alternatively, a depth can be associated with each pixel in each image, for example, by 
adding another dimension to each image and increasing the parameter set from 2D to 3D, 
and then, when combining images, the frontmost pixel 'from the images being combined is 
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the one that is used, as is done in 3D computer graphics, where a Z-buffer is used to select 
the frontmost pixels. Or alternatively, the exclusive or (XOR) of the first image and the 
second image can be used. 

S At step 412, a check is performed to determine whether the current output buffer is the last 
output buffer 218. If not, then the next output buffer, in this case the second output buffer 
216, is selected at step 414, and the process returns to step 404 to randomly select another 
M images from the V input buffers 208 to 212, and subsequently generate another 
superimposed output image 318 for storage in the second output buffer 216 by executing 

10 steps 406 to 410. These steps are repeated until all the V output buffers 214 to 218 have 
been used to store a new superimposed output image 318. Once the last output buffer 218 
has stored an image, the test at step 412 succeeds, and an iteration of the process is 
complete. Rather than selecting another output buffer, the contents of the output buffers 
214 to 218 are then used as input images to the next iteration of the process. This can be 

15 achieved by physically copying the contents of the output buffers 214 to 218 into the 
respective input buffers 208 to 212. However, it is more efficient simply to exchange, at 
step 416, the roles of the two sets of buffers by exchanging or swapping the corresponding 
buffer pointers into the random access memory 112. In any case, after the superimposed 
output images have been provided in respective input buffers 208 to 212 and the output 

20 buffers 214 to 218 have been cleared, the fractal generation process begins the next 
iteration by selecting the first output buffer 214 at step 402, and repeating the steps 
described above. As many iterations as desired can be performed. 

The contents of one or more of the output buffers 214 to 218 can be displayed on a display 
25 device associated with the fractal generation system for viewing by the user. Individual 
images selected by the user or images generated by a user-specified number of iterations of 
the fractal generation process can be stored as final output images 102 on the non-volatile 
storage medium 114 for subsequent use. After a number of iterations, the output images 
102 are a set of V-variable fractal images derived from the preceding input images. It will 
30 be appreciated by those skilled in the art that the although the image data generated by the 
fractal generation system is referred to as fractal image data, it may only represent a 
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nascent fractal if the finite number of iterations of the fractal generation process are not 
sufficient to produce a fractal. 

With successive iterations, the output images 318 become less dependent on the initial 
5 input images 104, and eventually become characteristic of the input parameter tables 106 
and their selection probabilities 108. The images 318 generated by the system then form a 
diverse sample of an infinite collection of all possible images associated with the initial 
parameter tables 106 and probabilities 108. In particular, because the transformations are 
contractive, the original input images 104 for the first iteration will eventually be 
10 contracted to the size of a single pixel when the entire output image is viewed, given a 
sufficient number of iterations. Thus the output images become independent of the 
original input images 104. Accordingly, the input images 104 can be alternatively 
provided or generated by the system itself, rather than being provided by the user as input, 
and the value of V can be provided as input to the system. By adjusting the parameter 
15 tables 106, the appearance and form of the classes of images generated by the system can 
be modified. For example, one class of images might resemble a group of tea trees, 
another a collection of faces, and another might resemble (and be used to represent) 
possible stock market simulations. As a result, the possible applications of the fractal 
generation system and process are very broad and diverse. 

20 

For example, Figures 6 to 10 are screenshots of images generated by the fractal generation 
system after 4, 5, 6, 1 1 and 44 iterations, respectively, using V=2 and the following input 
parameter tables: 



25 IFS 1: 





a 


b 


c 


d 


e 


f 


transformation 1 : 


0.5 


-0.375 


0.3125 


0.5 


0.375 


.1875 


transformation 2: 


0.5 


0.375 


0.1875 


-0.5 


0.375 


.6875 
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IFS 2: 





a 


b 


c 


d 


e 


f 


transformation 1: 


0.5 


-0.375 


0.3125 


-0.5 


-0.375 


.8125 


transformation 2: 


0.5 


0.375 


0.1875 


0.5 


-0.375 


.3125 



In this example, each transformation is an affine transformation defined by 
(x,y) -> (ax+by+c, dx+ey+f), where the parameters a to f are as indicated in the parameter 
5 tables. 

Figures 11 to 15 are screenshots of fractal images generated after 1, 3, 6, 8, and 22 
iterations, respectively, using V=2 and the following parameter tables: 

10 IFS1: 



a 


b 


c 


d 


e 


f 


e 


h 


I 


.151827 


.034074 


-.215212 


-.123706 


-.246197 


.054262 


-.212523 


.059302 


-.234997 


.011248 


.011321 


-.000425 


-.001814 


-.026793 


.062591 


.014512 


-.057923 


.086631 


S2: 

a 


b 


c 


d 


e 


f 


S 


h 


I 


.017380 


-.035039 


.046965 


.029178 


.006695 


.025663 


.021377 


-.012129 


.048469 


.144191 


.004070 


-.183496 


-.099432 


-.243409 


.043246 


-.151388 


.001065 


-.199968 



In this example, each transformation is a projective transformation defined by 
15 ( Xj y) _> ( (ax+by+c) / (gx+hy+I), (dx+ey+f) ) / (gx+hy+I) ) , where the parameters a to I 
are as indicated in the parameter tables. 

Figures 16 to 23 are screenshots of images generated after 1, 2, 3, 4, 6, 8, 22, and 40 
iterations, respectively, using the same projective transformation with V=2 and the 
20 following parameter tables: 
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IFS1: 



a 


b 


c 


d 


e 


f 


g 


h 


I 


.001231 


-.335970 


.002318 


-.164124 


-.028044 


-.125233 


.003548 


-.033609 


-.309261 


.002575 


.039210 


.007817 


-.045553 


-.005329 


.056137 


-.001851 


-.027965 


.126001 


S2: 






a 


b 


c 


d 


e 


f 


8 


h 


I 


-.006211 


-.197079 


.208860 


.087110 


.013677 


.011083 


-.008298 


.075200 


.205408 


.003412 


.306195 


-.303822 


.141253 


.023639 


-.286545 


.001890 


.028936 


-.312419 



5 The image data or images generated by the system represent a new type or class of fractals, 
referred to as V-variable fractals. As can be seen from Figures 6 to 23, they can be 
considered to be random fractals because they exhibit a degree of randomness, but this is 
intermediate between deterministic fractals and prior art random fractals generated by 
other methods. As such, they are able to represent natural forms more accurately than prior 

10 art fractals, and yet can be generated more readily than prior art random fractals. 

An image representing a V-variable fractal can be characterised as follows. Each fractal 
image is a combination of at least two smaller images. There is a set of V basic images 
such that each of the smaller images represents a contractive affine transformation of a 

15 basic image from the set of V basic images. The V basic images and the contractive affine 
transformations provide the first level of decomposition of the image representing the V- 
variable fractal. Each of the two or more smaller images from the first level of 
decomposition is itself a combination, of at least two smaller again images and there is a 
second set of V basic images such that each of the smaller again images represents a 

20 contractive affine transformation of a basic image from the second set of basic images. 
The second set of V basic images and the corresponding contractive affine transformations 
provide a second level of decomposition of the image representing the V-variable fractal. 
These decompositions can be continued for at least four levels of decomposition. However, 
if every set of basic images contains just one distinct basic image, then the representing 
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image is a deterministic or random fractal of a previously known type, and as such is not 
included in the new class of fractals. 

It will be apparent that a variety of alternative transformations can be used in the above 
5 manner. For example, projective transformations other than affine transformations can be 
used. 

Figures 6 to 10, 13 to 15, and 19 to 23 are examples of V-variable fractal images where 
V=2. Each of these images was generated by four or more iterations of the fractal 
10 generation process. At each scale there are at most V non-equivalent sets. These sets can 
be combined partially or wholly one on another. The sets typically change from 
magnification level to magnification level. While V=2 in Figures 6 to 23, larger values of 
V can be used. 

15 The properties of V-variable fractals are described in more detail in M. Barnsley, 
J. Hutchinson and 0. Stenflo, A New Random Iteration Algorithm and its Attractors, 
included as Appendix L 



20 The fractal generation system has been described above in terms of contractive affine 
geometrical transformation and superimposition of two-dimensional images. However, it 
will be apparent that the iterated function systems 106 can include other function systems 
that can be used to modify the input image data 104. For example, the iterated function 
systems 102 can include contractive projective transformations, or affine transformations 

25 that are contractive on average, where the average is determined by the probabilities 108 
associated with the iterated function systems 106. Because the system generates output 
images by combining multiple, transformed copies of input images, contractive 
transformations are preferably used in order to limit the output image size when translated 
and/or rotated images are combined, and this is preferably limited to the size of the input 

30 images. However, the transformations can alternatively be contractive on average, as 
described above. Similarly, it will be apparent that the transformed images can be 
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combined by methods other than superposition, as described above. Moreover, the fractal 
generation process can be applied to image data representing three-dimensional images to 
generate a new class of three-dimensional fractals, or indeed to input data of even higher 
dimensionality. The image data need not represent an existing image or be used to 
5 generate an image, but can be used for any desired purpose. 

In an alternative embodiment, the transformation modules 302, 304 can modify the colour 
and/or brightness of input images, in addition to the geometrical transformations described 
above. These additional modifications are defined by additional parameters included in the 

1 0 parameter tables 1 06. When the modifications of colour and/or brightness are performed in 
a contractive manner, the output images generated by the system, after a small number of 
iterations, depend only on the input parameter tables 106 and their selection 
probabilities 108. If the brightness of an image is represented by a number z, then an 
example of a contractive transformation of brightness is new_z= p*z + q where |p|<l and q 

15 is a constant. Contractive transformations of colour are achieved by three contractive 
transformations of brightness, one for each of the red, blue and green colour components. 

In a further embodiment, the parameter tables are generated graphically. For example, an 
affine transformation can be characterised by how it transforms a given triangle into 
20 another triangle. For example, a user can use a pointing device (such as a mouse) to drag 
the vertices of a given triangle to new positions and thereby define an affine 
transformation. 

In yet a further alternative embodiment, the parameter tables 106 and selection 
25 probabilities 108 are themselves generated from predetermined probability distributions, 
such as normal distributions, using a random number generator. 

It may be appreciated that the new class of V-variable fractals, realised as images or in any 
digital form, may also be generated by processes other than those described above. For 
30 example, those skilled in the art will realise that for the generation processes described 
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herein, which correspond to forward methods, there are also corresponding backward 
methods, although such methods will be less efficient for generating V-variable fractals. 

Many modifications will be apparent to those skilled in the art without departing from the 

5 scope of the present invention as herein described with reference to the accompanying 
drawings. 
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A NEW RANDOM ITERATION ALGORITHM AND ITS 
ATTRACTORS 

MICHAEL BARNSLEY, .JOHN HUTCHINSON, AND ORJAN STENFLO 

ABSTRACT. We describe a new class of random fractals and fractal measures 
with the remarkable property that they can be computed using a forward pro- 
cess; that is, a version of the Random Iteration Algorithm can be used to 
sample the random fractals. The results are obtained by constructing a new 
iterated function system (a super IFS) out of a collection of IFSs with prob- 
abilities. The attractor of the super IFS is a collection of sets or measures 
and associated (probability) measures, (what we call a superfractal.) When 
the underlying space is for example R 2 , and the transformations are compu- 
tationally straightforward, such as affine transformations, the superfractal can 
be sampled means of the algorithm, which is not only highly efficient in terms 
of memory usage but also the fractals that it produces are special; at each 
magnification level there is a predominance of a few key "forms" or "shapes" . 
The algorithm is illustrated by some computed examples. Some variants, spe- 
cial cases, generalizations of the framework, and potential applications are 
mentioned. 



1. Introduction and Notation 

1.1. Fractals and Random FVactals. A theory of deterministic fractal sets and 
measures, using a "backward'' algorithm, was developed in Hutchinson [5]. A dif- 
ferent approach using a "forward" algorithm was developed in Banisley and Deinko 
[12]. 

Falconer [10], Graf [17] and Mauldin and Williams [4] randomized each step in 
an approximation of a "backward" process to obtain random fractal sets. Arbeiter 
[21] introduced and studied random fractal measures; see also Olsen [23]. ** Insert 
a sentence on [7] and [8]. For further material see Zahle [25], Patzschke and Zahle 
[24], and the references in all of these. 

Deterministic fractal sets and measures are defined as the attractors of certain 
iterated function systems (IFSs), see Section 2 for details. Approximations in prac- 
tical situations quite easily can be computed using the associated random iteration 
algorithm. Random fractals are typically harder to compute because one has to 
first calculate lots of fine random detail at low levels, then one level at a time, build 
up the higher levels. 

In this paper we restrict the class of random fractals to ones that we call V- 
variable fractals. Superfractals are sets of V-variable fractals. They can be defined 
using a new type of IFS, in fact an IFS made of a finite number N of IFSs, and there 
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is available a novel random iteration algorithm: each iteration produces a new V- 
variable fractal, lying increasingly close to the superfractal, and moving ergodicaliy 
around it. 

Superfractals appear to be a new class of geometrical object, their elements lying 
somewhere between deterministic fractals which correspond to V = N = 1, and 
random fractals which correspond to V = oo. They seem, to allow geometric mod- 
elling of some natural objects, examples including reaiistic-iooking leaves, clouds, 
and textures; and good approximations can be computed fast in elementary com- 
putergraphical examples. They are fascinating to watch, one after another, on a 
computer screen, diverse, yet ordered enough to suggest coherent natural phenom- 
ena and potential applications. 

Areas of potential applications include computer graphics, stock market simula- 
tion, and simulation using finite V of a case with V = oo corresponding to Brownian 
motion. 

1.2. An Example. Here we give an illustration of an application of the theory in 
this paper. By means of this example we introduce informally V-variable fractals 
and superfractals. We also explain why we think these objects are of special interest 
and deserve attention. We start with two pairs of contractive affine transformations, 
IfJ./J} and {/?,/?}, where /£ : □ □ with □ := [0,1] x [0,1] C R 2 . We use 
two pairs of screens, where each screen corresponds to a copy of □ and represents 
for example a computer monitor. We designate one pair of screens to be the Input 
Screens, denoted by (Di.Cb). The other pair of screens is designated to be the 
Output Screens, denoted by (Oi'tCk')' 

Initialize by placing an image on each of the Input Screens, as illustrated in 
Figure 2, and clearing both of the Output Screens. We construct an image on each 
of the two Output Screens as follows. 

(i) Pick randomly one of the pairs of functions {Si , /<&} or {/? ,/|}, say^j 11 }. 
Apply /" l to one of the images on Di or U 2j selected randomly, to make an image 
on Dir. Then apply /J 1 to one of the images on Oi or D 2 , also selected randomly, 
and overlay the resulting image J on the image now already on Di. (For example, 
if black-and-white images are used, simply take the union of the black region of J 
with the black region on Oi, and put the result back onto U\.) 

(ii) Again pick randomly one of the pairs of functions {/},/}} or 8av 
{/r 3 i/? 2 }- Apply /" 2 to one of the images on U u or Q 2 , selected randomly, to 
make an image on Dy. Also apply f$ 7 to one of the images on n u or Cfe, also 
selected randomly, and overlay the resulting image on the image now already on 
□2/. 

(iii) Switch Input and Output, clear the new Output Screens, and repeat steps (i), 
(ii) and (iii) many times, to allow the system to settle into some sort of "stationary 
state". 

What sorts of images do we see on the successive pairs of screens, and what sorts 
of images are produced in the "stationary state"? What sorts of things does the 
theory developed in this paper tell us about such situations? As a specific example, 
let us choose 

(1.1) fx (*. V) = (0.5x - 0.375y + 0.3125, 0.5s + 0.375y + 0.1875), 



(1.2) /t (*»V) = (°- 52 + °' 375 2' + °- 1875 ' -°' 5x + °' 375y + °' 6875 )' 
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Figure 1. Triangles used to define the four transformations 
AS/2 1 , /?, ««d /?• 

(1.3) fi(x,y) = (0.5* - 0.375y + 0.3125, -0.5* - 0.375y + 0.8125), 

(1.4) fi(x, y) = (0.5x + 0.375y + 0.1875, 0.5x - 0.375y + 0.3125). 

We describe how these transformations act on the triangle ABC in the diamond 
ABCD, where A = (0.25,0.5), B = (0.5,0.75), C = (0.75, 0.5), and D = (0.5,0.25). 
Let Bi = (0.28125,0.71875), B 2 = (0.71875,0.71875), B 3 = (0.28125,0.28125), 
and B 4 « (0.71875,0.28125). See Figure 1 (to be inserted). Then we have 

(1.5) ft (A) = A, fl{B) = B U fl(C) = B; 

(1.6) /i(A) - B, &{B) = B 2 , ${C) = C; 

(1.7) /?(A)=A^(B) = B3,/?(C) = B; 

(1.8) /2 2 (A)=A/?(B) = B 4 ,/|(C) = C. 

In Figure 2 we show an intial pair of images, two jumping fish, one on each of the 
two screens Oi and Cfe. In the subsequent Figures, 3, 4, 5, 6, 7, 8, and 9, we show 
the start of the sequence of pairs of images actually obtained in a particular trial, for 
the first seven iterations. Then in Figures 10, 11, and 12, we show three successive 
pairs of screens, obtained after more than twenty iterations. These latter images 
are typical of the ones obtained for any number more than twenty iterations, very 
diverse, but always representing continuous "random" paths in K 2 j they correspond 
to the "stationary state", at the resolution of the images. Notice how the two 
images in Figure 11 consists of the union of shrunken copies of the images in Figure 
10, while the curves in Figure 12 are made from two shrunken copies of the curves 
in Figure 11. 
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Figure 2. An initial image of a jumping fish on each of the two 
screens di and Cb. 




FIGURE 3. The pair of images after one iteration. 




Figure 4. The two images after two iterations. 
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FIGURE 5. The two images after three iterations. Both images are 
the same. 





Figure 6. The two images after four iterations. Both images are 
again the same, a braid of fish. 





FIGURE 7. The two images after five iterations. The two images 
are the same. 



MICHAEL BARNSLBY, JOHN HUTCHINSON, AND 6RJAN STBNPLO 




FIGURE 10. Images on the two screens Di and after a certain 
number L > 20 of iterations. Such pictures are typical of the 
"stationary state" at the printed resolution. 
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FIGURE 11. Images on the two screens after L + 1 iterations. 




Figure 12. Images on the two screens after L -h 2 iterations. 



This example illustrates some typical features of the theory in this paper, (i) 
Keeping attention on the sequence of images that occur in one screen, we see a 
sequence of images, one per iteration, (ii) After sufficient iterations for the system 
to have settled into its " stationary state", each image looks like a finite resolution 
rendering of a (different) fractal set, each fractal belonging to the same family, in 
the present case a family of continuous curves.(iii) In fact, it follows from the the- 
ory that the pictures in this example correspond to curves with fractal dimension 
between one and two, with this property: for any e > 0 the curve is the union of 
"little" curves (i.e ones such that the distance apart of any two points is no more 
than e) each of which is an affine transformation of one of at most two continuous 
closed paths in K 2 . (iv)The successive images, or rather the abstract objects they 
represent, are related as follows: all lie upon a single object, itself a component of 
an IFS attractor of a certain IFS in a certain space. This component is called a 
superfractal, the images are distributed upon it according to the invariant proba- 
bility measure of the IFS, and the images produced by the algorithm are distributed 
according to this measure, (v) The images produced in the "stationary state" are 



8 



MICHAEL BARNSLEY, JOHN HUTCHINSON, AND 6RJAN STENPLO 



independent of the starting images. For example, if the initial images in the exam- 
ple had been of a dot or a line instead of a fish, and the same sequence of random 
choices had been made, then the images produced in Figures 10, 11, and 12 would 
have been the same at the printed resolution. 

1.3. The structure of this paper. The main contents of this paper, while con- 
ceptually not very difficult, involves potentially elaborate notation because we deai 
with iterated function systems (IFSs) made of IFSs and probability measures on 
spaces of probability measures. So a material part of our effort has been towards 
a simplified notation. Thus, below, we set out some notation and conventions that 
we use throughout. 

The core machinery that we use is basic IFS theory, as described in [5] and 
[12]. So in Section 2 we review relevant parts of this theory, using notation and 
organization that extends to and simplifies later material. To keep the structural 
ideas clear, we restrict attention to IFSs with strictly contractive transformations 
and constant probabilities. Of particular relevance to this paper, we explain what 
is meant by the random iteration algorithm. 

(** Check later.) We begin Section 3 with a discussion of the code trees that 
underlie random fractals. We show that there exists a special class of trees, called 
V-variable trees, where V is an integer. Unlike general trees, these objects can be 
characterized in terms of an IFS acting on an appropriate space. What are these 
trees like? At each level they have at most V distinct subtrees! The IFS enables 
one to put measure attractors on the set of such trees, such that the trees can be 
sampled by means of the random iteration algorithm. 

(♦♦Check later.) The discussion of trees in the first few subsections of Section 
3 is recapitulated several times over in the later subsections, where the same basic 
IFS theory is applied in more and more elaborate settings. We arrive at the formal 
concepts of V-variable fractals and super-fractals. 

** In Section 4 we describe some generalizations and potential applications of 
the theory. 

1.4. Some Notation. We use notation and terminology consistent with [12]. 
Throughout we reserve the symbols M, N, and V for positive integers. We will 

use the variables m e {1,2, ...,M}, n € {1,2, ...,7V}, and v € {1,2,...,V}. 

Throughout we use an underlying metric space (X,dx) where the space is X and 
the metric is d — dx* Except where otherwise stated, (X,dx) is assumed to be a 
compact metric space. We write % v to denote the compact metric space 

(1.9) X x X x ... x X. 

V TIMES 

with metric 

(1.10) d(x,y) = d x v(x t y) = Max {d x (x v ,y v ) : v = 1,2,..-, V}, Vx,y 6 X M , 

where x = (xi,X2» .~xv) and y = (yi,l/2i >»yv)* 

In some applications, to computer graphics for example, (X, d\) is a bounded 
region in R 2 with the Euclidean metric. We will be concerned with continuous 
transformations / : X X, for example affine and projective transformations acting 
on R 2 . 
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Let S = §(X) denote the set of all subsets of X, and let C € S. We extend the 
definition of a function /:X->Xto/:S-»Sby 

(1.11) f{C) = {/(*) 

Let H=1HI(X) denote the set of non-empty compact subsets of X. Then if / : X — > 
X is continuous we have / :H -» EL We use dm to denote the HausdorfF distance on 
M. implied by the metric d% on X. This is defined as follows. Let A and B be any 
pair of sets in H, define the distance from A to B to be 

(1.12) V(A>B) = Mox{Min{d x (*,y) : V € B} : x € A}, 
and define 

(MS) d H {A,B) - Max{V{A, £), V(B, A)}. 

Then (M,d H ) is a compact metric space. We will write (M v ,dwr) to denote the V- 
dimensional product space constructed from (IIMa) just as (X v , d^y) is constructed 
from (X,dx). When we refer to continuous, Lipshitz, or strictly contractive func- 
tions acting on M v we assume that the underlying metric is d H v - 

We will in a number of places start from a function acting on a space, and extend 
its definition to make it act on other spaces, while leaving the symbol unchanged 
as above. 

Let B a B(X) denote the set of Borel subsets of X. Let P = P(X) or M = M(X) 
(depending on the intended interpretation as indicated below) denote the set of 
Borel probability measures on X. In some applications to computer imaging one 
sets X = [0, 1] x [0, 1) C R 2 and identifies a black and white image with a member of 
IHI(X). Greyscale images are identified with members of M(X). Probability measures 
on images are identified with P(H(X)) or P(M(X)). 

Let djp (X ) denote the Monge-Kantorovitch metric on P(X). This is defined as 
follows. Let \i and u be any pair of measures in P. Then 

d P ( Ml i/) = Sup{ I /d/i- / fdu : / : X - K, \f(x) -/(y)| < d x (x,y) V x,y € X}. 
X X 

Then (P,dp) is a compact metric space, as is (M,d M ). We define the push-forward 
map/:P(X)-P(X) by 

(1.15) f(ji) = M o r x V /* € P(X). 

Again here we have extended the domain of action of the function / : X — » X. 

We will use such spaces as P(M V ), and P((M(X)) K ) (or M(1HI V ), and M{M V ) 
depending on the context). These spaces may at first seem somewhat Baroque, but 
as we shall see, they are very natural. In each case we assume that the metric of 
a space is deduced from the space from which it is built, as above, down to the 
metric on the lowest space X, and often we drop the subscript on the metric without 
ambiguity. So for example, we will write 

(1.16) d(A, B) m dbcm>" )< A > B ) V A » 3 6 M « M ( X )) V )- 
We also use the following common symbols: 

N = {1,2,3, ...}, and X = {... - 2, -1, 0, 1, 2, ...}. 
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2. Iterated Function Systems 

2.1. Iterated Function Systems. In this section we review relevant information 
about IFS. To clarify the essential ideas we consider the case where all mappings 
are contractive, but indicate in later sections how these ideas can be generalized. 
The machinery and ideas introduced here are applied repeatedly later on in more 
elaborate settings. 
Let 

(2.1) F= {X;/i,/ 2 ,...,/m;Pi,P2. -.Pm} 

denote an IFS with probabilities. The functions f m : X -> X are contraction map- 
pings with fixed constant 0 < I < 1; that is 

(2.2) d(/m(*), /m(y)) < I ■ d(x,y) Vx,y € X, Vm € {1, 2, M}. 

The p m 's are probabilities, with 

M 

(2.3) £>m = 1, Pm > 0 Vm. 

We define mappings F : IHI(X) EI(X) and F : P(X) — P(X) by 

M 

(2.4) F(K) = (J f m {K) VJT 6 H, 

m=l 

and 

Af 

(2.5) F(M) = £ Pm/mW V/X € P. 

m=l 

In the latter case note that the weighted sum of probability measures is again a 
probability measure (greyscale image). 

Theorem 1. [5]The mappings F : H(X) M(X) and F : P(X) -> P(X) are both 
contractions with factor 0 < / < 1. That is, 

(2.6) d(F(K), F(L)) < I • d(K, L)VK,Le IH(X), 
and 

(2.7) rf(F(Ai),F(i/)) < I • t/) V i/ € P(X). 

.As a consequence, there exists a unique nonempty compact set A 6 H(X) such that 

(2.8) F(A) = j4, 
and a unique measure p. 6 P(X) such that 

(2.9) F(M) = /*. 

7%e support of fi is contained in, or equal to A, with equality when all of the 
probabilities p m are strictly positive. 

Definition 1. The set A in Theorem 1 is called the set attractor of the IFS F, and 
the measure ft is called the measure attractor of F. 
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We will use the term attractor of an IFS to mean either the set attractor or the 
measure attractor. We will also refer informally to the set attractor of an IFS as 
a fractal and to its measure attractor as a fractal measure. Furthermore, we say 
that the set attractor of an IFS is a deterministic fractal. This is in distinction to 
random fractals, and K-variable random fractals which are the main goal of this 
paper. 

There are two main types of algorithm for the practical computation of extrac- 
tors of IFS that we term deterministic algorithms and random iteration algorithms, 
also known as backward and forward algorithms, c.f. [15]. These terms should not 
be confused with the type of fractal that is computed by means of the algorithm. 
Both deterministic and random iteration algorithms may be used to compute deter- 
ministic fractals, and as we discuss later, a similar remark applies to our V-variable 
fractals. 

Deterministic algorithms are based on the following: 
Corollary 1- Let A 0 e Bfl(X), or^e P(X), and define recursively 

(2.10) A k = F^fc-O, or » k = F(ji M ) 9 Vfc € N, 
respectively; thenlS 

(2.11) LimAk = A, or Lirnfi k = fi, 
respectively. The rate of convergence is geometrical; for example, 

(2.12) d(A ky A)<l k - d{Ao, A) V* G N. 

Example 1, An example of a deterministic algorithm involving a simple IFS on 
the unit square in R 2 . See Figure IS. 

The random iteration algorithm is based on the following theorem, see [10]. 

Theorem 2. Specify a starting point x\ 6 X. Define a random orbit of the IFS to 
be {x t }fl x where xj +1 = f m (x t ) with probability p m • Then for almost oil random 
orbits for any xi € X, we have: 

(2.13) M(B)= wK£ni£i^i£iH } . 

for all B e B(X) such that p(BB) = 0, where BB denotes the boundary of OB. 

Remark 1. This is equivalent by standard arguments to the following: for almost 
all random orbits the sequence of point measures ±(5 Xl + <5* 2 + ... + S Xt ) converges 
to fx, see for example [13], pages 11 and 12. 

Example 2. An example of the random iteration algorithm involving a simple IFS 
on the unit square in R 2 . See Figure U 

The following theorem expresses the ergodicity of the IFS F. The proof depends 
centrally on the uniqueness of the measure attractor, and it applies more generally, 
as stated for example in [11] (??) where a probabilistic proof is given. 

Theorem 3. Suppose that \l is the measure attractor for the IFS F and B e B(X) 
is such that f m (B) C B Vm € {1, 2, M). Then a{B) = 0 or 1. 
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FIGURE 13. An illustration of the deterministic algorithm, also 
showing the texture effect. 




Figure 14. "Picture" of the measure attractor of an IFS with 
probabilities produced by the random iteration algorithm. The IFS 
is {□;//, /2»^2i0.4, 0.3,0.4} where the transformations are those 
defined in the Introduction. The measure is depicted in shades of 
green, from 0 (black) to 255 (bright green). 
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Proof. Let us define the measure p\B restricted to B) by (fi[B)(E) = fi(BC\E). 
The main point of the proof is to show that p[B is invariant under the IFS F . (A 
similar result applies to /i|_B c where B° denotes the complement of B. 
If E C B c , for any m, since f m (B) C B, 

(2.14) fMB)(E) » MB n /->(£)) - M«) = o. 

Moreover, for any m, 

(2.15) KB) = fMB)(X) = /m(MlB)(B). 
It follows that 

M 

KB) = ^PmfmKB) 



m=l 

Af M 



(2.16) = jP Pm/m(MlB)(B) + £ PnJm(M|S C )(S) 

msl m=l 
M 

= MB) + J! Pm/m(ML^ C )(^) (from (2.15)). 

T7l=l 

Hence 

M 

(2.17) £pm/«(MlB c )(fi)=0. 

Hence for any measureable set E C X 
(H[B)(E) - Ai(B0£7) 

M 

(2.18) = £> m / m /i(BnB) 

mr=l 

Af M 

- £ p m fMB){B OE)+J2 PmfMB c ){B n f?) 

M 

= Ep^^^)(^)+0 (using (2.16)). 



Thus /i[B is invariant and so is either the zero measure or for some constant 
c > 1 we have cfi[B = p (by uniqueness)^ fi[B + p[B c . This implies fi[B c = 0 
and in particular A*(B C ) = 0 and /x(B) = 1. □ 

Theorem 3 provides a simple explanation for an computer graphics effect that 
may be seen when the set attractor of an IFS on R 2 is computed using a deter- 
ministic algorithm, as illustrated in the following example. A similar effect can also 
occur in the case of V-variable fractals. 

Example 3. Example illustrating the texture effect 

**An important number concerning fractals is the fractal dimension. Not only 
is it an important invariant, but in practice says something about the look and feel 
of an attractor. 

Definition 2. **Fractal Dimension 
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Definition 3. ** Open set condition 
Theorem 4. **Fractal Dimension 

2.2. Code Space. A good way of looking at an IFS F as in (2.1) is in terms of 
the associated code space E = {1,2,...,M}°°. Members of E are infinite sequences 
from the alphabet {1,2, M} and indexed by N. We equip E with the metric d £ 
defined for u x by 

(2.19) d £ (w,>f) = i 

where k is the index of the first symbol at which w£ x differ. Then (E,di:) is a 
compact metric space. 

Theorem 5. Let A denote the set attractor of the IFS F. Then there exists a 
continuous onto mapping F : E — > A, defined for all 0x0*20*2— € E by 

(2.20) F(ai<r2<T3-) = 0 f°* 0 - 0 (*)• 

The limit is independent ofxeX and the convergence is uniform in x. 

Definition 4. The point am**... € E is called an address of the point F(a\OQP&~) € 

A. 

Note that F : E — *• A is not in general one-to-one. 

The following theorem characterizes the measure attractor of the IFS F as the 
push-foward, under F : E A t of an elementary measure p € IP(E), the measure 
attractor of a fundamental IFS on £. 

Theorem 6. For each m € {1, 2, M) define the shift operator s m : E -+ E by 

(2.21) SmfoWa. ») = nur^^a- 

V<ri<72^3— € E. Then s m is a contraction mapping with contractivity factor \. 
Consequently 

(2.22) S:= {E;s 1| 52>—i«m;P1iP2i— iPAf} 

is an IFS. Its set attractor is E. Its measure attractor is the unique measure p € 
P(E) such that 

(2.23) = Poi • P02 ' ••• " P^». 

V* € N, Vcr lv cri,... f <7* € {1,2,..., M}. 

If pis the measure attractor of the IFS F, with F : E— ► A defined as m Theorem 
5, then 

(2.24) H = F(p). 

We call 5 the shift IFS on code space. It has been well studied in the context 
of information theory and dynamical systems, see for example [14], and results can 
often be lifted to the IFS F. **Example. 
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FIGURE 15. Illustration of a 2-variable tree, up to level 5. 



3. Trees of Iterated Function Systems and SuperFractals 

3.1. TVees of IFSs. Let (X,dx) be a compact metric space, and let M and N be 
a positive integers. For n 6 {1,2, TV} let F" denote the IFS 

(3.1) {X; ft> P?, ...pXf } 

where each /£ : X X is a Lipshitz function with Lipshitz constant 0 < I < 1 and 
the p£i 's are probabilities with 

M 

(3.2) ^ P ; = l,p;>OVm,r l . 

m=l 

Let 

(3.3) r = {X; F\ F 2 , F"; P 1 , F 2 , .... P N }> 

where the F n, s are probabilities with 

N 

(3.4) ^F Tl = l,F n >0 Vn. 

F< n > > 0 Vn € {1, 2, TV} and ^ (r ° = We ^ ^ a for masons 

that will become clear as we proceed. This is not a standard IFS because the 
constituent "functions", the F n 's, themselves standard IFS's, are to be composed 
in a special way, to be explained. 

We associate with T the following collection of labelled trees ft. Let T = Tm 
denote the (M-fold) tree of finite sequences from {1,2, ...,M}, including the empty 
sequence 0. A 2-fold tree is illustrated in Figure 15. For cr = a\o^a k 6 T 
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let |a| = k. If also r = TiTv^.ti e T then err is the concatenated sequence 
cr 1 <r2...<yfcrir2...ri. A code tree is a map a : T -> {1, 2, -AT} and we write 

(3.5) H - {cr|<r : T - {1, 2, JV}} 

for the set of all such code trees. 

We define a metric on fl by, for w ^ x, 

(3.6) = ~ 

where A; 'is the index of the first symbol at which the strings that represent w and 
>c disagree. Then (£2,cfo) is a compact metric space. 

A construction tree of T is a code tree wherein the symbols 1,2, and N are 
replaced by the respective IFSs F l ,J ?2 , and F N . A construction tree consists 
of nodes and branches, where each node is labelled by one of the IFSs belonging 
to We will sometimes associate the M branches that lie below a node with the 
constituent functions of the IFS of that node; taken in order. 

We use the notation to denote the set of construction trees of T* For a 6 ft 
we write to denote the corresponding construction tree. We will use the same 
notation T{o) to denote the random fractal set associated with the construction 
tree .F(<7), as described in the next section. 
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3.2. Random Fractals. In this section we describe the canonical random fractal 
sets and measures associated with the superlFS T in (3.3). These sets and measures 
are defined using random code trees obtained as follows. We put a probability 
measure p e P(ft) on the set of code trees, by setting <x(r) = n with probability 
P n > for n € {1,2,..., N] independently at each node of the tree r. € T. Thus we 
are able to speak of the random code trees ft with probability distribution p $ and of 
selecting trees a € ft according to p. 

Let a superlFS T be given as in (3.3), let k € N, and define 

(3.7) F h : ft x H(X) -+ H(X), 
by 

(3.8) Fk(<r){K)= |J /# 0 > o J*** o ... o f*r***-"-*(K) 

V a 6 ft and K € H(X). (The set F k (p)(K) is obtained by taking the union of 
the compositions of the functions occuring on the branches of the construction tree 
F(c) starting at the top and working down to the fc th level, acting upon K.) 

In a similar way, with measures in place of sets, and unions of sets replaced by 
sums of measures weighed by probabilities, we define 

(3.9) ^:ftxP(X)->P(X) 
by 

(3.10) ^(a)(r,)= Y, ° # T ° o ... o 

{reT M \\T\=k} 

V a e ft and r/ € P(X) where 

(3.11) = 

V m € {1, 2, M} and n 6 {1, 2, iV}. Note that the JF fc (a)(T/) all have unit mass 
because of the in Equation 3.11. 

Theorem 7. Let sequences of functions{F k } and {fk} be defined as above. Then 
both the limits 

(3.12) H°) = Lim{F k (<r)(K)}, and #(<r) « Zim{^ fc (cr)(7;)}, 

exist, independently ofK and v, and the convergence (in the Hausdorff and Monge- 
Kantorovitch metrics, respectively Jis uniform in a, K, andn,. The resulting func- 
tions 

(3.13) r : ft-> H(X) and f : ft-*- P(X) 
are continuous* 

Proof Make repeated use of the fact that, for fixed a € ft, both mappings are 
compositions of contraction mappings of contractivity l t by Theorem 1. □ 

" Let 

(3.14) fy = {^(a) € H(X)|cr € ft}, and h - {^) €E P(X)|a € ft}. 
Similarly let 

(3.15) V = ^(rf « p o F~ x € P(fi), and*- F{p) ^po^e P(ij). 
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Definition 5. The seta S) and h ore called the sets of random fractal sets and 
random fractal measures, respectively, associated with the superJFS T. These ran- 
dom fractal sets and measures are said to be distributed according to ip and Vp, 
respectively. 

Random fractal sets and measures are hard to compute. There does not appear 
to be a general simple forwards (random iteration) algorithm for practical compu- 
tation of approximations to them in two-dimensions with affine maps, for example. 
The reason for this difficulty lies with the inconvenient manner in which the shift 
operator acts on trees a € Q relative to the expressions (3.8) and (3.10). 

Theorem 8, ** Fractal dimension of random fractals, say in the case of simili- 
tudes. 
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3.3. Contraction Mappings on Code Trees and the Space (lv V € N. 
This parameter will describe the variability of the trees and fractals that we are 
going to introduce. Let 

(3.16) n v = ft x Q x ... x a 

V TIMES 

We refer to an element of fi v as a grove. In this section we describe a certain IFS 
on QY t and discuss its set attractor Sly: its points are (V-tuples of) code trees that 
we will call V -groves. 

One reason that we are interested in fly is that, as we shall see later, the set of 
trees that occur in its components, called V -trees, provides the appropriate code 
space for a V-variable super facta!. 

We are going to use here and in later sections an index set 

(3.17) A := {1, 2, N} v x {1, 2, .... V} MV . 
A typical index a € A will be denoted 

(3.18) a = {nijiv n v ) x (vi,»2>.«i»Mv) 

where each 6 {1, 2, N) and each Vj € {1, 2, V}. 

We are also going to need a set of probabilities {i^la € ,4}, with 

(3.19) ^2 P° = 1, P a > 0 Va € A 

a€-4 

These probabilities may be more or less complicated. Some of our results are 
specifically restricted to the case 

(3.20) P* = . 

which uses only the set of probabilities {Pl l \P&\ ...,P (/V) } belonging to the su- 
perlFS (3.3). 

Theorem 9. For each n 6 {1,2, N} define the n th shift mapping ( n : Q M ft 
by 

(3.21) (£„(w))(0) = n and £ rt ("))(mr) = w m (r) V r € T, m € {1, 2, M}, 

foru- (wi,w 2 , u/Af ) € H M . ZViae is, fte mapping f n creates a corfe free with its 
top node labelled n attached directly to the M trees wi,W2i ..., w Af • For each a€A 
define Va := Q v ->Q V by 

(3.22) ?n 2 (Wv M+ i « w « M+a) -i w va M )» - 

T/ien ij a : fi v — ♦ £2 V is a contraction mapping with Lipshitz constant £ and conse- 
quently 

. (3.23) $ := {Q v \ rj a9 P° ,c€ .4} 

w an /F5 0/ strictly contractive maps. 

Let the set attractor and the measure attractor of the IFS $ be denoted by IV, 
and vl v respectively. Clearly, Qy € M(Q V ) while fiy e 1?(Q V ). The elements of toy 
are certain V-tuples of M-fold code trees on an alphabet of N symbols, which we 
characterize in Theorem 12. But we think of them as groves of trees. 
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Definition 6. Sly is called the space of V-groves. An element of Sly « called a 
V -grove. An element of a component of a V -grove is called a V-tree. 

For all v € {1,2, V), let us define Sly tV C Sly to be the set of v th components 
of groves in Sly. 

Theorem 10. For allv€ {1,2, V) we have 

(3.24) Qy tV = *V,i := {set of all V-trees). 

When the probabilities {P°\a e A} obey (3.20), the distribution of trees w € SI that 
occur in the v ih components of groves produced by the random iteration algorithm 
corresponding to the IFS $ converge weakly to the marginal probability measure 

(3.25) fi Vil (B) := n v (B, SI, SI, Q)VB €= B(ft), 
independently of v. 

Example 4. Describe how this works in practice. Give an example, with one code 
twice as probable as the other. 

Theorem 11* The set of V-trees converges to the set of all trees SI in the metric 
ofM(Sl),thatis 

(3.26) £wnJV.i = ft. 
Moreover, when the probabilities {P°|o € A] obey (3.20), 

(3.27) = P, 

where p is the stationary measure on trees introduced in Section 3.2, and conver- 
gence is in the metric ofH(Sl). 

Proof Let V > M k . Then both 

(3.28) d^a^Sly.uSl)^^, 
and 

(3.29) d v{{l) (u Vil ,p) < ^y. 

□ 

Theorem 12. Let a € Sl v . Then a € Sly if and only if the number of distinct 
subtrees of the trees in the V -grove cr, at any fixed depth, is at most V. Given 
any tree a e SI with the property that, at any depth, it, contains at most V distinct 
subtrees, then there is a V -grove cr 6 ft y whose first component is a. 

Let 5V = A 00 * This is simply the code space corresponding to the IFS From 
Theorem 5 there exists a continuous onto mapping $ : £y Sly defined by 

(3.30) $(aia 2 a 3 ...) = ° Va 2 0 - ° la* ( w ) 

for all aia 2 a 3 ... € for any w S ft. 

In the terminology of section 2.2 the sequence 010203- € £v is an address of the 
V-grove #(a l a203...) € Sly and £y is the code space for the set of V-groves Sly. 
However, as an address space for indexing the V-groves of Sly it is not particularly 
helpful because in general $ : £v -» Sly is not one-to-one as various compositions 
of its functions can coincide, as the following example shows. 
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Example 5. Example showing that $ : Ey -> &v is not 1-1 in general. 
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3.4. Contraction Mappings on M v and the SuperfVactal Set Sy.i- 
Theorem 13, Let V 6 N, let A be the index set introduced in (3.11), let T be 
given as in (3.3), and let probabilities {P a \a € -4} be given as in (3.19). Define 

(3.31) ?*:M V -+M V 
by 

(3.32) f*{K) = ( (J /£(tf v J, (J &(K VM + m ) 9 ... (J C(^v«.. + J) 

V/ir=(iirx,i^2,....^v)€lffi v ,Va€A 2%en 

(3.33) = {H v ;^°,P a ,a € .4} 
is an IFS with, contractivity factor I. 

Proof. We begin by noting that, V{K U K 2 ,„.,K M ), (L u L<i,...,L M ) e M M , 

M[J ft(K m ),\J Kit***)) 

(3.34) < Max{d^fZ(K m ),f^(L m ))} 
< Max{l • dn(K m , L m )} 

= l-das,(K,L). 
Hence, V(K lt K2,...,K v ),(L u L2,...,Lv) <=M V , 

<V (^(tfi, tf 2 K v ), ^°(Li, La, .... L„)) 

(3.35) = Afoa:{dH( U fZ(K VltM _ M+tn ), (J /£•(£„„„_„+„,))} 

m=sl mssl 

< A40X{Z •^ HW ((^« V A/-/M+l»^W«M-Ai+a»"*»^ l> M)» 

< < ■ d HV ((Ku .... Kv), (^i ■ I*, 

□ 

The theory of IFS in section 2.1 applies to the EPS Ty. It possesses a unique set 
attractor 9)y € M(H V ), and a unique measure attractor tyy € P(M ). The ran- 
dom iteration algorithm corresponding to the IFS Ty may be used to approximate 
sequences of points (vectors of compact sets) in Sjy distributed according to the 
probability measure tyy. 

However, the individual components of these vectors in %y, certain special sub- 
sets of X, are the objects we are interested in. Accordingly, for all v € {1, 2, V}, 
let us define %y tV C EE to be the set of v th components of points in fiy. 

Theorem 14. For all v € {1, 2, V} we have 

(3.36) S}y tV = flm. 

When the probabilities in the superlFS Ty are given by (3.20), the sets K € H that 
occur in the V th component of vectors produced by the random iteration algorithm 
converge weakly to the marginal probability measure 

(3.37) VvAB) := Vv(B v H t H > ..^B)VB 6 B(H), 



RANDOM IFSS 



23 



independently of v. 

Definition 7. We caUS)v t i a superfnxctal set Points in S)v,i are called V-variable 
fractal sets. 

Example 6. An example that will be referred to below. Points, each one a fractal 
set appear to dance around on the superfractal. 

By Theorem 5 there is a continuous mapping Tv - %v — ► that assigns to 
each address in the code space Sv a V-tuple of compact sets in f)v But this 
mappings is not helpful for characterizing S)y because Tv ' £v ~ * &v 13 not in 
general one-to-one, for the same reason that # : Ev — ♦ Hy is not one-to-one, as 
demonstrated in Example 5. Also Ev does not exhibit in any direct manner the 
tree-like structures that underlie V-variable fractals. The following is closer to the 
point. It tells us that the set of K-trees is a useful code space for V-variable fractals. 

Theorem 15. Let Sl v ,v denote the set ofv th components of members ofSlv Let 
T : H— » 1HI(X) be the mapping introduced in Theorem 7. Then 

(3.38) T(ttv,i)=Si VA . 

When the probabilities in the superlFS Tv are given by (8.20) we have 

(3.39) T(ji Vtl )=VPvA 

where fi V i marginal probabUity distribution given by (3.25). 

Definition 8. The code tree ^(ayo^a^...) is called a tree address of the V-variable 
fractal T{aia 2 a3...). 

The mapping T : fiy.i -» fiv.i together with Theorem 12 allows us to charac- 
terize V-variable fractals as follows: 

Theorem 16. Any V-variable fractal is made of at most V "forms" or "shapes" 
at each level 

Example 7. Example of the Theorem, in the form of a black and white picture, 
pointing out the key fact! 

Theorem 17. The set of V-variable fractal sets associated with the superlFS T 
converges to the set of random fractal sets introduced in Section 3.2; that is, in the 
metric of M(M(X)), 

(3.40) ^i™ 55 ^ 1 = 

Moreover, when the probabilities {P°|a 6 A} obey (3.20), then in the metric of 
P(M(X)) 

(3.41) ^jnV^i = 

where is the stationary measure on random fractal sets associated with the su- 
perlFS T 

Proof. This follows at once from Theorems 7 and 11. □ 
Theorem 18. Result about fractal dimension, say in the case of open set condition. 



24 



MICHAEL BARNSLEY, JOHN HUTCHINSON, AND OR J AN STENFLO 



3.5. Contraction Mappings on M v and the Superfractal Measures S) Vt u 
Let M v » (M(X)) V (= (fP(X)) v '). In this section we follow the same lines as in 
Section 3.4, constructing an IFS using the individual IFSs of the superEPS 
except that here-the underlying space is M v instead of H^. 

Theorem 19. Let V 6 N, let A be the index set introduced in (3.17), let T be 
given as in (3.3), and let probabilities {i^jfi € .4} be given as in (3.191 Define 

(3.42) F*:M V -+M V 

by 

(3.43) 

M M M 

- (£ pmoo, £ icr/ff&w)."-- £ c/s v ^v-v + J) 

M= (MhA*2i-iMv)^ mV - Tnen «F° : M v M v is a contraction mapping and 

(3.44) = {M v ;^° ,.P a ,a € .4} 
i5 an IFS with contractivity factor I. 

The set attractor of the IFS Ty is hv e H(M V ), a subset of M v , a set of 
V-tuples of probability measures on X. As we will see, each of these measures is 
supported on a V-variable fractal set belonging to the superfactal 3v,i- Tne measure 
attractor of the IFS fy is a probability measure % v € P(M V ), namely a probability 
measure on a set of V-tuples of normalized measures, each one a random fractal 
measure. The random iteration algorithm corresponding to the IFS Fy may be 
used to approximate sequences of points in hv> namely vectors of measures on X, 
distributed according to the probability measure 

As in Section 3.4, we define hv t v to be the set of v th components of sets in %y, 
for v€ {1,2, ...,K}. 

Theorem 20. For all v € {1, 2, V} we have 

(3.45) hv>v = hv t \. 

When the probabilities in the IFS Fy are given by (3.20), the probability measures 
p € rVH(X) that occur in the V th component of points produced by the random iteration 
algorithm converge weakly to the marginal probability measure 

(3.46) *v,i(B) := *u(B,M,M,...,M)VB € B(M) f 
independently of v. 

Definition 9. We call hv t x * superfractal set of measures (of variability V). Points 
in hv % \ ore called V-variable fractal measures. 

Example 8. An example that will be referred to below. Points, each one a fractal 
measure, appear to dance around on the superfractal set of measures. 

By Theorem 5 there is a continuous mapping Fy : IV j"> hv that assigns to 
each address in the code space Ey a V-tuple of measures in hv- But this mappings 
is not helpful for characterizing hv because Fy : Sv — » hv is not in general one- 
to-one, for the same reason that $ : Ev -> SV is not one-to-one, as exhibited in 
Example 5. 
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Theorem 21. Let Slv t v denote the set of V th components of members of Sly Let 
F : S7-» IP(X) be the mapping introduced in Theorem ??. Then 

(3.47) ^W = flv,i« 

When the probabilities in the superlFS T are given by (3.20) we have 

(3.48) Hnv % i) = %v 9 \ 

where (jl vi is the marginal probability distribution given by ($.25). 

Definition 10. The code tree ^(a^a^...) is called a tree address of the V-variable 

fractal measure ^V( a i a 2 a 3--)- 

The mapping F : SV f i -» %v t i together with Theorem 12 allows us to charac- 
terize V-variable fractals as follows: 

Theorem 22. Any V-variable fractal measure is made of at most V "forms" or 
"shapes" at each level 

Example 9. Example of the Theorem, in the form of a greyscale picture, pointing 
out the key fact! 

Theorem 23. The set of V-variable fractal measures associated with the superlFS 
T converges to the set of random fractal measures introduced in Section 3.2; that 
is, in the metric o/IEfl(M(X)) 

(3.49) , Ifmflm = h- 
Moreover, when the probabilities {P°Ia e A} obey (3.20), 

(3.50) £i™^ 1== *' 

where the convergence is in the metric o/P(M(X)), and where % is the stationary 
measure on random fractal measures associated with the superlFS T . 

Proof This follows at once from Theorems 7 and 11. D 
Theorem 24. ** Fractal dimension theorem, say in the case of similitudes. 



26 



MICHAEL BARNSLBY, JOHN HUTCHINSON, AND ORJAN STENFLO 



4. Algorithms, Applications and Examples 
Here we consider some of the theory and applications.. 
4.1. Applications. 
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